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ABSTRACT 

We present a sample of 17 newly discovered ultracool dwarf candidates later than ^M8, drawn from 
231.90 arcmin 2 of Hubble Space Telescope Wide Field Camera 3 infrared imaging. By comparing the 
observed number counts for 17.5 < J125 < 25.5 AB mag to an exponential disk model, we estimate a 
vertical scale height of z sc l = 290 ± 25 (random) ±31 (systematic) pc for a binarity fraction of ft, = 0. 
While our estimate is roughly consistent with published results, we suggest that the differences can 
be attributed to sample properties, with the present sample containing far more substellar objects 
than previous work. We predict the object counts should peak at J125 ~ 24 AB mag due to the 
exponentially-declining number density at the edge of the disc. We conclude by arguing that trend 
in scale height with spectral type may breakdown for brown dwarfs since they do not settle onto the 
main sequence. 

Subject headings: Keywords: Galaxy: structure — Galaxy: stellar content — stars: low-mass, brown 
dwarfs 



1. INTRODUCTION 

Star counts have long been used to determine the 
structure of our Galaxy. Early attempts were plagued 
by patchy extinction and by mathema tical instab ilities 
in the inversion of star counts (see iBol H93l . but 
IBahcall fe Soneiral (|1980l I1984D revived the endeavor 
by avoiding regions with significant extinction, and 
by fitting a physically motivated model with only a 
few parameters (see Bahcall 1986 for a review). Re- 
cently, there has been renewed interest in star counts as 
ever-cooler stella r populations are discovered and need 
to be modeled (IRvan et all 120051: iPirzkal et al.l " 
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Pirzkal et al. 2009; IDeacon et aHl2009tlBochanski et all 
2010t iDelorme et all 120101) . Because brown dwarfs cool 
and change spectral type on relatively short timescales, 
their vertical scale heights may reflect not just Galactic 
structure, but also their cooling times. 

The Galactic distribution of the ultracool dwarf pop- 
ulation has garnered much interest from a community 
studying far more distant objects. Since the extremely 
red optical and near-infrared colors of the ultracool 
dwarfs are similar to tho se of Lyman-break galaxies 
(LBGs; ISteidel et al.l Il99fif ) at redshifts of 5 < z < 7, 
there are concerns that the high-redshift galaxy sam- 
ples may be contaminated by these Gal actic objects (e.g. 
iCaballeroT Burgass er. fc Ele ment 2008) . In the absence 
of spectroscopic identification, LBG studies often resort 
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to statistically correcting t heir number counts fo r fore- 
ground contamination (e.g. iBouwens et all 120061 ). Nat- 
urally, this correction relies on accurately characteriz- 
ing the G alactic distr i bution of ultracool dwarfs. To 
this end, Ry an et ail ([2005D identify 28 dwarfs with 
{%' — z')ab > 1.3 mag (which are types later than ~M6; 
Bochanski et al. 2010) in 15 parallel fields from the Ad- 
vanced Camera for Surveys (ACS) aboard the Hubble 
Space Telescope (HST). B y assuming an exponential disk 
model. [Rvan et al~l (|2005l ) derive a vertical scale height of 
z sc i = 350 ± 50 pc, and conclude that the de epest surveys 
of zc^ Q LBGs were >97% pure. Similarly, Pirzk aTet all 
(|2005f) find a vertical scale height of z scl = 400 ± 100 pc 
from three spectroscopically identified lat e-M and early- 
L dw arfs in the Hubble Ultra Deep Field (jBeckwith et all 
120011 . 

The overwhelming majority of ultracool dwarfs to 
date have been found in shal l ow, very wide-field sur- 
veys Ce.g. iDelfosse et all 119991: iKirkpatrick et al.l 119991: 



2004) , and more recently with deepe r 



iKnapp et al.l 

datasets (e.g. IDelorme et all 12001 IDeacon etall 120091 ). 
Since t hese objects are int rinsically very faint (M,' > 
16 mag: lHawlev et all 120021 nearly all know n ultracool 
dwar fs reside within ~100 pc of the Sun (e.g. lReid et al.1 
2008), which makes determining the Galactic-scale dis- 
tribution difficult or impossible. While this issue can 
be mitigated to a large extent by probing further into 
the disk, limited observing time and detection efficiency 
have restricted studies t o narrow fields-qf-view and/or 
single lines -of-sight (e.g. IRvan et al.l 120051 : IPirzkal et all 
120051 12009D . Naturally this leads to simplified models, 
large uncertainties on model parameters, and significant 
variations between authors. In this paper, we begin to 
overcome these limitations by drawing our sample from 
very deep parallel and pointed fields with HST, which 
have the sensitivity to find an LO-dwarf to ~3.5 kpc and 
a TO-dwarf to ~ 700 pc. These represent a significant 
increase in survey distances for ultracool dwarfs. 
This work is organized as follows: in §[2]we describe the 
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observations and source catalogs, in § [3] we discuss our 
ultracool dwarf sample selection, in § [4] we present our 
analysis and scale height measurement, in §[5] we assess 
our systematic uncertainties, and in §[6] we conclude with 
a brief review and thoughts toward future improvements. 
Unless explicitly stated otherwise, all magnitudes an d 
colors are given in the AB system (|Oke fc Gundll~983f ). 

2. OBSERVATIONS 

Here we discuss the parallel and pointed fields with 
HST which consititute our dataset. We list their salient 
properties in Table [TJ 

2.1. The HST Parallel Imaging 

The bulk of the data analyzed here come from the 
high-level science products from the Hubble Infrared 
Pure Parallel Im aging Extragalactic Survey (HIPPIES3; 
lYan et all 120111) . At present, this survey combines two 
pure parallel imaging programs with HST Wide Field 
Camera 3 (WFC3; PropIDs: 11700 and 11702) and coor- 
dinated parallels^ from the Cosmic Origins Spectrograph 
(COS) guaranteed time observations (GTO). Every field 
has infrared imaging in F098M, F125W, and F160W, 
and optical imaging in F606W and/or F6OOLF0. By de- 
sign, these fields are at relatively high Galactic latitude 
(|6| > 20°), have a total exposure time of >4 HST orbits, 
and sample random pointings through the Galaxy. 

The data reduction and mo saicking of the HIPPIES 
data are discussed in detail bv lYan et al.l (|2011[ ). but we 
will mention the key steps relevant for this work. Stan- 
dard procedures were followed except for enhanced re- 
moval of image defects. Since the HST parallel data 
are rarely dithered, the affected pixels were corrected 
by interpolating over neighboring pixels with the FIXPIX 
routine in IRAF. The main side-effect of this procedure 
is to decrease the usable area of each parallel point- 
ing. In total, we analyzed 39 parallel fields which cover 
185.41 arcmin 2 . Finally, we note that the COS GTO 
parallel fields also have additional imaging in F300X, 
F475W, and F475X, though we place no constraints on 
the colors in these bands. 

2.2. The WFC3 Early Release Science Program 

In addition to the parallel data, we include the WFC3 
Early Release Science program (ERS; PropID: 11359) 
data taken in the southern field of the Great Observa- 
tories Origins Deep Survey (GOODS-S; iGiavalisco et al.1 
l200l . The WFC3 imaging in the F098M, F125W, and 
F160W-bands augments the existing optical data from 
the ACS and covers 46.49 arcmin 2 . The ERS data have 
at least double the exposure time of the WFC3 paral- 
lel data in all bands. Details of the ERS data, such as 
experimental de sign, reduction, and im aging properties 
are discussed bv lWindhorst et al.l (|2011l ). Like the COS 
GTO parallels, the ERS subset of the GOODS-S field has 
been observed in many additional optical and ultraviolet 
bands, however we will not impose constraints on those 
colors to ensure a uniformly selected sample. 

7 http:/ /archive. stsci.edu/prepds/hippies/ 

8 See the HST User Information Report UIR-2008-001 for a dis- 
cussion of parallel imaging with HST. 

9 The fields at 02 h 20 m , 07 h 50 m , and 12 h 09 m have both optical 
bands. 
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Fig. 1.— HST/WFC3 filter throughput curves. The five filters 
used here are shown as F606W (V606; solid blue), F600LP (^600; 
dotted blue), F098M (Y 09s ; solid cyan), F125W (J125; solid green), 
and F160W (-H160; solid red). In light gray lines, we show select 
spectra from the A. Burgasser SpeX compilation (see Tabled with 
infrared types indicated on the left. The COS coordinated parallels 
and the ERS field have additional blue filters from HST. However, 
we do not explicitly place constraints on these colors, in order to 
ensure a uniform sample. 

2.3. Photometry 

We measu re all magnitude s as MAG_AUT0 with 
SExtractor (jBertin fc Arnoutsl 119961 ) in dual -image 
mode using the F125W image for detection. We require 
a minimum area of 5 connected pixels, a threshold (per 
pixel) of > 1.5 o~ over the local background, and use a 
5x5 pix Gaussian filter with full width at half maximum 
of 2 pix for source detection. All images are drizzled 
to / /090 pix -1 . We use 64 deblending sub-thresholds, 
a minimum contrast parameter of 10~ 4 , and a clean- 
ing efficiency of 10. We adopt the AB zeropoints from 
iKalirai et all (|2009aHH ) for the WFC3 data of F606W = 
26.08 mag, F600LP = 25.85 mag, F098M = 25.68 mag, 
F125W = 26.27 mag, and F160W = 25.96 mag, and 
F606W = 26.486 mag for the ACS data in the GOODS- 
S field. Based on these object catalogs, we compute the 
50% completeness depth in the F125W image as the mag- 
nitude where the ratio of the observed counts to a power- 
law fit falls to 0.5, and report these depths in Table Q] 

The parallel and ERS data give us a total of 40 in- 
dependent sight lines through the disk which cover a 
total of 231.90 arcmin 2 . As mentioned above, each 
field has the same three near-infrared bands; however, 
the optical imaging differs in wavelength and instru- 
ment. In Figure [1] we show our optical and infrared 
bandpasses: F606W (solid blue), F600LP (dashed blue), 
F098M (cyan), F125W (green), and F160W (red), which 
we refer to as V606, ^600, ^098, ^125, and #160, respec- 
tively. 

3. ULTRACOOL DWARF CANDIDATES 

3.1. Sample Selection 

To ensure that our objects are point-like, we require 
the axis ratio to be (b/a) > 0.7 and the half-light ra- 
dius as measured by FLUX_RADIUS in SExtractor^ to 

10 SExtractor will measure the radius at which a some fraction 
of the total flux is reached based on the setting PHOT.FLUXFRAC, 
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TABLE 1 
Survey Fields 



r 1G1Q 


D At 

rh m s\ 
{ ) 


Ttor>T 
JJCC 1 

{ ) 


■t 1 

[aeg) 


fit 



(dcg) 


AO? 

l u J 


T* 

J 5Q 

(mag) 


A** 
A J 

(mag) 


wpticai 
joana 


parUHU— UzZZ 


m 1 n nn a k 
Ul 1U Uy.4o 


no oo oo n 
— UZ ZZ Zo.U 


1 Q O (1Q70')0 

loo.yo/ ZOZ 


ft/1 C/IOICO 

— d4.o4Z1oz 


4.68 


27.69 


0.04 


rSU 1 H 


^ .-no 1 q i i o k /i 

parUzlo+lzo4 


no 1 Q QO 7C 

UZ lo OO.(0 


1 1 O C A KO O 

+iz o4 oy.z 


1 KO fl1 07/10 

lOZ.Ulo 1 4Z 


A K Oft 1 Q Qft 

— 4o.zbloob 


4.69 


27.02 


0.11 


TTftnnT TD 

r bUULF 


COSUZZ / — 41U1 


no 07 Kft m 

UZ 2 / ou.yi 


A 1 m Q A a 

—41 Ul o4.4 


O K A 1 ft 1 Oftn 

Zo4.1bloby 


ftK 7fl070n 

—bo. /yz / zy 


4.69 


28.00 


0.01 


TTftnnT to 
r bUULiF 


_ _ _ n o a n 10C7 
cosUZ4U— loo / 


no /in 07 ft o 
UZ 4U Z / .bo 


1 O K T 1 A A 

— loot 14.4 


onn ft/i nnnn 

zuu.b4yuuy 


ftO ftQftVQn 

—bo. bob /oU 


4.69 


27.95 


0.03 


TTftnnT to 
r bUULiF 


ERS 


no o o o o no 

Uo oz Zo.Uo 


TT A O K n O 

— Z7 42 5U.Z 


ooo a r\'vnrzr\ 

ZZo. 4U7959 


K/i /i/ii/ino 

— o4.4414Uo 


46.49 


28.20 


0.01 


1 bUbW 


cosU4oy— oolb 


n/i Qfi occ /io 
U4 oy Z0.4Z 


co -\ /? ACv A 

— oo lb 4U.4 


Oft 1 Q o A n A o 

Zbl.oo4y4o 


/in n a ft OT'ft 
— 4U.y4bZ /b 


4.69 


28.21 


0.00 


TTftnnT "D 
r bUULiF 


_„„nc;Q O ft/ino 

paruooy— b4uy 


nc: on on CO 

uo oy ou.oz 


ft A no OQ /l 
— b4 Uy Uo.4 


07Q ftKO'7/1 r 7 

Z / o.boU /4 / 


QO m K/l 7fl 

— oZ. Ulo4 /U 


A 70 

4. /Z 


Oft A K 

Zb.40 


n nK 
U.Uo 


T?ftnft \S7" 
r bUb W 


parUooo— b4Uo 


nc ro nft no 
Uo oo Ub.UZ 


ft /i ok i o n 
— b4 Uo lo.U 


070 KOKftftO 

z / o.ozobbo 


Qn KOKKKT' 

— oU.ooooo / 


/I 7ft 

4. /b 


27.01 


0.04 


TTftnft'XST' 

1 bUbW 


parUbzo— b4ol 


nft o q q /i nft 
Uo zo o4.Ub 


a a o 1 /i n 1 
— b4 ol 4y.i 


ov/i ooonn/i 
Z / 4.ZoZyy4 


OT* Oft A O A ft 

— Z / .Zb4Z4b 


4.68 


26.32 


0.05 


r bUbW 


parUbzo— o4v>y 


nc OQ /IO 1 Q 

Ub zo 45. lo 


ft a on a i n 
— b4 oy 41. U 


07 '/l QOOftOT* 

ZY4.oozbo / 


OT* O K Q 70n 

— z /.zoo /oU 


/I 71 

4. / 1 


26.94 


0.05 


T7 , ftnftl!(7' 

r bUbW 


„ _ ..nft q 7 7tri n 

parUbo r— /oiy 


nft 07 fir no 

Ub o / Uo.Uz 


TK 1 o on A 

— to lo oy.4 


OQft /i 1 nnnn 
zob.4iyuUU 


07 (170 1 ft 1 

— Z / .U /olbl 


6.92 


26.98 


0.09 


T7 , ftnftl!(7' 

r bUbW 


parU/oU+zyi/ 


n7 Kn Kn kq 

U/ oU oU.oo 


i on ift k o ft 
+zy lb oo.b 


101 QKQOQ/1 

iyi.ooooo4 


I O/l (lCfl9A7 

+Z4.ybUoU/ 


4.81 


27.43 


0.04 


d(J 1 H 


paru / 0o+oU4o 


n7 kk c'7 no 
U / OO / .Uo 


1 Qfl A Q 1 n o 

+oU 4o iu.y 


1 On 01 A QOft 

iyu.zi4oyb 


1 Oft A KQ KO*7 

+zb.40ooy / 


A ftQ 

4. bo 


Oft 71 

Zb. / 1 


n nft 
U.Ub 


T?ftnft WT 

r bUb W 


parUoUo+oy40 


no no 01 qq 
Uo Uo Zl.oo 


I QO A K OK Q 

+oy 40 ZO.o 


1 QO OOQ K/1 /I 

iou.yzoo44 


I Q 1 1 oo Km 
+01.1Z00U1 


A ftQ 

4. bo 


OK A O 

zo.4y 


n c\a 
U.U4 


T?ftnft WT 

r bUb W 


paruoiy+4y 11 


no m in n/i 
Uo iy iy.U4 


i /in 1 1 OK A 

H~4y ii uo.4 


1 7fl OOQftVQ 

i /u.uyob / o 


1 QA O/l A QQft 
+ 04.Z44oob 


A ftQ 

4. bo 


07 70 

2( . t 2 


n nK 
U.Uo 


T?ftnft \S7" 
r bUb W 


noon i nooo 

parUozU+zooz 


no on no a 1 

Uo zU Uo. 41 


i oo o o n k n 

-f Zo oz Uo.U 


"i nn ooonKO 

lyy.ozoyoz 


1 on oomno 

- i r 2\).o2b2>}2 


4.69 


27.29 


0.04 


1 bUbW 


,-.-^,-,nQ /I ft 1 7ft K '3 

cosUo4b+ /boo 


no /1ft OO 

Uo 4b ZZ.ob 


1 'zft KO on o 

+ to oo oy.o 


1 oft ftn'T'n'7'7 

lob.bU/y/ / 


1 Q O 7Cm Q K 

+oz. /bUloo 


4.68 


28.35 


0.02 


TTftnnT to 
r bUULiF 


„ _ „nnnc; i no k k 
parUyUo+Uzoo 


nn n k q? k o 
uy Uo o / .OZ 


i no K K 1 ft 

+Uz oo ol.b 


OOft Q A Q 1 7Q 

ZZb.o4ol / o 


i on nft 1 n a 
+oU.yblU4o 


4.68 


27.03 


0.03 


r bUbW 


paruyuy— uuui 


no no no i a 

uy uy uy.i4 


nn m /i*7 1 
— UU Ul 4 / . 1 


OQO Q1COQ1 

ZoU.oloUol 


1 QO 1 O/l OQft 

+ou. iy4yob 


A ftQ 

4. bo 


07 Qft 

2 1 .ob 


n oq 
U.Uo 


T?ftnft WT 

r bUb W 


_._nni /i i oqoi 
paruy 14+ZoZl 


no i/i ift qo 

uy 14 ib. oz 


1 OQ 01 A A a 
+Zo Zl 44. b 


1 OQ 1 A VOOft 

iyo. 14 / uzb 


1 /in QC KOQ /I 

+4Z.OOOyo4 


/I ftQ 

4. bo 


07 KO 

2 1 .0U 


n no 
U.UZ 


T?ftnft WT 

r bUb W 


«nY.noQi i /iKnK 
paruyzi+40Uo 


no 01 3B 1 K 

uy zi oo.io 


i a k nK oq n 
+40 Uo Uo.U 


1 VK 1 A OQ1 n 

1 / 0.14Z01U 


1 /i /i (inm on 
+44.yUUlZU 


A ftQ 

4. bo 


07 1 ft 

z / . lb 


n no 
U.UZ 


T?ftnft WT 

r bUb W 


_ _ ,,nno k i a A ok 
paruyzo+44Zo 


nn ok q o ik 

uy zo oz. io 


I A A O K 1 O 

+44 zo ol.o 


i /o.yoy4y4 


1 /IK ft /I Q 1 'zn 

+4o.b4oi /y 


4.69 


27.85 


0.01 


TTftnnT "D 
r bUULiF 


paruyzo+4UUU 


no OK Q K /IK 

uy zo 00.40 


i a n nn 1 q n 
+4U UU lo.U 


1 CO QOI OftO 

loZ.oZlUbZ 


I /IK Q'T'QQQft 

+40. /ooob 


A ftQ 

4. bo 


07 A O 

z / .4y 


n ni 
U.U1 


T?ftnft \S7" 
r bUb W 


parlUoU+ooUo 


in on ko ko 
1U oU OZ.OZ 


i q o no O A K 
+00 Uo Z4.5 


1 no KftKQO^ 

loo.obooz / 


1 K Q ftft K 1 1 

+oo.bbolzl 


4.68 


27.62 


0.01 


Trftnc^sr 
1 bUbW 


cosllol+oll f 


11 q 1 on no 

ii oi zy.yo 


I Q 1 17 OI O 
+ol It Zl.O 


1 n a 7'}onoQ 
iy4. /ozUoo 


i 70 nn /1007 

+ /z.uy4oz / 


4.69 


27.98 


0.02 


TTftnnT to 
r bUULiF 


_ _ _ 1 1 ri i K /I /1 1 

pari lo 1+0441 


11 k 1 /in oft 

ii oi 4y.zb 


I K A A n KO O 

+54 4U oy.o 


1/in /iQKnQ/i 
14U.4ooyo4 


1 ftn Q 70ft OK 

+bU.o / zbzo 


A 71 

4. / 1 


27.70 


0.01 


r bUbW 


~, _ „ i -i pro i nnKft 
pal HOZ+UUob 


1 1 KO A Q OO 

ii oz 4o.yz 


i nn kk ki o 
+UU 00 01. Z 


070 OOCOI Q 

Z / Z.ZZoZlo 


1 ftO OK K1 OK 

+bU.Z0olZ0 


/I ftQ 

4. bo 


07 7Q 

2( . to 


n no 
U.UZ 


T?ftnft \~\t 
r bUb W 


-i onn i /ic/in 

parlzUy+454o 


i o nn o a o o 

Iz Uy z4.oz 


1 /IK A o or 1 

+45 4d ZD.l 


AAA O C7C C C 

144.3676D6 


+69.615667 


4.70 


28.09 


0.01 


BOTH 


parlz44+ooob 


1 O A A /IK oi 

lz 44 45. Zl 


I O O K c n E T 


134.455667 


1 OO n j4 O A A 1 

+oo.U4o441 


4.68 


28.13 


0.01 


h 606 W 


oni nnnn 
parloUl— UUUU 


iom ift ft 1 
lo Ul Ib.bl 


nn nn it n 
— UU UU z r .U 


3U8. 312235 


I ft O 7ft 1 O A 7 

+bz. /blo4/ 


4.68 


27.38 


0.02 


r 0UUL,F 


parl336— 0027 


13 36 48.75 


— OU 27 57.9 


326.341678 


+60.326874 


4.68 


27.86 


0.03 


F600LP 


parl340+4123 


13 40 31.87 


+41 23 U3.3 


90.813506 


+72.543346 


4.68 


28.17 


0.01 


F600LP 


parl436+5043 


14 36 56.72 


+5U 42 58.6 


89.753150 


+59.068438 


4.69 


28.19 


0.01 


F606W 


parl524+0954 


15 24 10.17 


+09 54 19.8 


14.751734 


+50.137021 


4.68 


27.64 


0.04 


F600LP 


parl631+3736 


16 31 34.28 


+37 36 21.4 


60.300146 


+43.026120 


4.68 


27.83 


0.01 


F606W 


parl632+3733 


16 32 18.38 


+37 33 24.3 


60.246932 


+42.877428 


4.68 


27.64 


0.01 


F606W 


cos2057-4412 


20 57 22.01 


-44 12 26.9 


356.582832 


-40.624781 


4.68 


27.45 


0.03 


F600LP 


cos2202+1851 


22 02 48.68 


+18 5U 58.5 


76.653695 


-28.493144 


4.92 


28.10 


0.06 


F600LP 


par2345-0054 


23 45 02.34 


-00 54 11.0 


88.894697 


-59.313985 


4.68 


27.83 


0.03 


F600LP 


cos2350-4331 


23 50 36.39 


-43 31 30.3 


335.844576 


-69.509649 


4.68 


28.18 


0.01 


F600LP 



t Coordinates refer to the field center in the J2000 epoch. 

t Solid angle in arcmin 2 . 

*The approximate 50% completeness limit. 

"The J-band extinction from Sclilcgcl, Finkbcincr, & Davis (1998). 



be 1.2 < rso < 1.8 pix. While unsaturated point- 
sources should have half-light radii roughly independent 
of brightness, we find a weak trend in the stellar locus, 
therefore our exact half-light radii limits vary slightly 
with magnitude. Using the morphological criteria pre- 
sented in Figure [2J we identify 5982 point sources in the 
231.90 arcmin 2 surveyed. 

We select our ultracool dwarf candidates from the cat- 
alog of point sources based on their optical and near- 
infrared colors. We synthesize empirical (V(306 — ^098), 
(ieoo- 5^098)) (*098 -«/x25)j and (^125 — #160 ) colors from 
a library of spectra from the 3-meter NASA Infrared 
Telescope Facility compiled by A. BurgasseiEB (listed 
in Table 0) . In Figure [31 we show the infrared color- 
color diagram for known L-dwarfs (green triangles), T- 



dwarfs (red circles), and M-dwarfs, giants, and subdwarfs 
(blue, magenta, and cyan squares, respectively). Un- 
fortunately, our broadband data cannot accurately con- 
strain the spectral type since the J125- and iii6o-bands 
equally sample the strong H 2 absorption at A = 1.34 /mi 
(see Figure [T]). Therefore, the only unique spectral types 
we can derive from these near-infrared data are given by 

MLT f 0.0 < (J125 - Hi &0 ) < 0.65 mag; and , , 
\0.75 < (Yogs - J125) < 1.0 mag 1 > 

t f 0.0 < (J125 - ffi6o) < 0.65 mag; and 
1 1.0 < (F 098 - J 160 ) < 0.7 x (J 125 - ff 160 ) + 1.0 maV 



which we adopt as 05; r , , , 

11 Distributed at http://web.mit.edu/ajb/www/browndwarfs/spexprism/ T I - 125 _ Ul 60) - U ' U ma S' anC1 

LJLU — " ' 1 1 0.75 < (Y 125 - J 160 ) < 1.0 mag 



(3) 
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TABLE 2 
SpeX Catalog" 1 " 



Reference 


Number of 




Citations 



A. Burgasser (unpublished) 
Burgasser et al. (2010 ) 
Burgasser et al. (2004a) 
Chiu et al. (20061 
Burgasser^et^l^ (2008a) 
Burgasser^e^il^ (2006a) 
Muench et al. (20071 

LooDcr. Kirkpatrick, & Burgasser (2007a) 

Burgasser. Burrows. & Kirkpatrick (2006b) 

Burgasser (2007a ') 

Siegler et al. (2007) 

Bur gasser et al. (2007b 1 

Burgasser fc McElwain (2006) 

Sheppard fc Cushing (20091 

Cruz et al. (20041 

Burgasser & Kirkpatrick (2006) 

Looper et al. (20081 

Burgasser et al. (2004b) 

Looper et al. (2007b 1 

McElwain fc Burgasser (20061 

Reid et al. (20061 

Muno et al. (20061 

Burgasser et al. (2008b1 

Burgasser et al. (20091 

Luhman et al. (20071 

(2006) 



Kirkpatrick et al. 
Burgasser (2007c 1 
Burgasser (2007d1 
Licbcrt fc Burgasser (2007) 



389 
116 
87 
51 
45 
19 
17 
11 
12 
8 
7 



* Compiled A. Burgasser and distributed at 

http: / / web. mit.edu /aj b/ www /browndwarfs /spexprism / 




20 22 
J 126 (AB mag) 

Fig. 2. — Morphological selection criteria. We use the half-light 
radius, axis ratio, and brightness as measured by SExtractor to 
identify point sources. The small dots represent all objects from 
the 231.90 arcmin 2 analyzed here, and the solid lines show the 
stellar locus selection region. We find a very weak relationship 
between half-light radius and brightness. Our constraints on half- 
light radius and axis ratio are consistent with known stars selected 
from the Sloan Digital Sky Survey. 
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Fig. 3. — Color-color diagram for ultracool dwarf selection. We 
show library of empirical L- and T-dwarfs from the SpeX spectro- 
graph compiled by A. Burgasser (see Table [2} as green triangles 
and red circles, respectively. We show the M-dwarfs, giants, and 
subdwarfs as blue, magenta, and cyan squares, respectively. Based 
on these colors, we define unique spectral types as equations JTJ, 
Ig), and (3). 
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Fig. 4. — Color-color diagram for ultracool dwarf candidates. 
We show the selection region defined by equations |[TJ, ©, and J3} 
as thick lines and all objects passing our brightness and morpho- 
logical criteria as small points. In blue, green, and red hatches, 
we show the expected colors of M-, L-, and T-dwarfs to highlight 
the ambiguity in spectral typing these stars with these bandpasses. 
The points with uncertainties represent our 17 ultracool dwarf can- 
didates in Table [3] 

For all types, we require (V606 — Yogs) > 2.0 mag and/or 
(^600 — Yq9s) >0.5 mag. We present our ultracool dwarf 
candidates in Table [3l 

3.2. Contaminants 

In addition to the ultracool dwarfs, there are three ad- 
ditional types of known astrophysical objects which may 
satisfy our optical and near-infrared color criteria: early- 
type galaxies (ETGs), LBGs, and high redshift quasars. 
In the absence of spectroscopic confirmation, we can only 
make statistica l arguments on these potential obje cts: 

Based on the lColeman. Wu. &: Weedman (1980) spec- 
tral templates, we estimate that an ETG at 1.2 <z<l.5 
will have optical and infrared colors similar to our ultra- 
cool dwarfs. By extrapolating the luminosity function 
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of ETGs at z = 1 (jFaber et all 120071 ). we find that our 
fields could have ~ 80 ETGs at these redshifts. How- 
ever, the number of ETG contaminants in our sample 
should be much lower, since we require our stellar can- 
didates be unresolved, which is not represented in this 
simple brightness calculation. To estimate the fraction of 
these ETGs that are also unresolved, we perform a simple 
Monte Carlo simulation. We draw 10 5 random absolute 
magnitudes from the z — 1 ETG luminosity function over 
our sample brightness range, which we convert to stellar 
masses assuming the mass-to-light of = 1 M Q Lq . 
For a given stellar mass, we draw a random size according 
to the mass-size relation for local ETGs from the SDSS 
(|Shen et al.ll2003l ). and determine the measured effective 
radius by quadratically adding the size of the Ji25-band 
PSF r mc 







SDSS 



Finally, we take the fraction 



of deviates which satisify our size criterion (see Figure [5]) 
as the fraction of detectable ETGs which would be unre- 
solved in these HST images. We estimate the poten tial 
ETG fraction to be <0.1% for the lShen et all ((20031) re- 
lation, and this fractio n only rises to 5, 1 % for the z^2 
mass-size relations (e.g. Ryan et al. 2010). Therefore, we 
conclude that our sample is largely free of contaminating 
ETGs. 

Like the ETGs, LBGs and quasars can only corrupt our 
sample in a very specific redshift range of 6.8<z<7.2. 
While both populations are likely to be unresolved, LBGs 
are typically far too faint, and quasars are far too rare, 
to have been included in our sample. In our bright- 
ness range (17.5 < J125 < 25.5 mag), we expect to 
find ~ 0.02 LBGs and ~ 0.01 quasars assu ming the 
iBouwens et"aH (|2010f ) and iWillott et all (|2010D luminos- 
ity functions, respectively. Therefore, we conclude our 
sample is likely free of any LBGs and/or quasars. 

4. ANALYSIS 

We determine the vertical scale height of ultracool 
dwarfs by comparing our observed star counts to those 
predicted from a Galactic structure model. We model 
the spatial distribution of dwarfs as 



n(r,z)oce- (r - r ° )/r - ! e- |z|/ ^ 



(4) 



where r Q = 8 kpc is the Solar position , r sc ; = 2 kpc is the 
radial scale length (|Juric et al.l I2008T) . and the constant 
of proportionality is set by the local (R < 20 pc) lumi- 
nosity funct ion, <fr(M). We ta ke the empirical luminosity 
fu nctions oflCruz et aLl (|2007t ) for the M8-L9-dwarfs, and 
of iRevle et alj (|2010l ) for T-dwarfs. The model number 
counts for the ith field are 

p 00 

N m (£ h bi)dm = ACliCi(m)dm / R 2 n(r l )^(M) dR, 



(5) 

where (ti,bi) are the Galactic coordinates, Af2j is 
the solid angle subtended, M = m — 51og(i?) — 
5 — A(£i,bi,R) is the absolute magnitude, R is 
the heliocentric distance in parsecs, A(£i,bi,R) is 
the extinction (discussed in more detail in § 15. ip . 



5 + R? cos 2 bi — 2Rtq cos cos 6, is the distance 
along the Galactic midplane, Zi— Rshxbi is the distance 
above the Galactic midplane, rj = \J x\ + zf is the Galac- 
tocentric distance (Bahcall 1986), and Ci(m) is the com- 
pleteness. The total model number counts are given by 
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Fig. 5. — The completeness corrections for the parallel (solid 
gray) and ERS field (dotted black). We estimate these corrections 
as the recovery rate of 1000 randomly-placed point sources in each 
magnitude interval. The completeness corrections are generally 
> 90% for our magnitude range and rarely tend to 100% (even at 
the bright-end), due to source blending. 



the sum over all of the fields 



N m = 



iVfielcU 

E 

i=i 



N m (li,bi), 



(6) 



which is parameterized by the vertical scale height in the 
Galactic model. We estimate the magnitude-dependent 
completeness by placing point-sources of known bright- 
ness at random locations within each field, cataloging the 
images as discussed in § 12.31 and taking the complete- 
ness as the fraction of recovered objects. In this way, 
we encapsulate the effects of our choice of SExtractor 
settings and source blending. In Figure [5j we show the 
completeness corrections for each parallel field (light gray 
lines) and the ERS field (dashed black line). The hatched 
regions indicate the magnitude ranges that we omit in 
our analysis. Given our relatively conservative limits of 
17.5 < J125 < 25.5 mag, these completeness values are 
rarely < 90%, and generally do not fall to their half- 
maximum values until Ji25~26 mag. 

We compute the model number counts for 240 < z sc i < 
350 pc with Az sc i = 5 pc and maximize the likelihood 
of obtaining the measured counts. Since the observed 
counts are in the limit of small integers, the probability 
distribution should be modeled as a Poisson distribution 
(|Cashlll979D . To avoid confusion with the usual Gaussian 
probability distributions, we denote this goodness-of-fit 
statistic as C 2 , and maximize the likelihood (£) in the 
usual way: 



C 2 (z sd ) = -21n 




N„ 



(7) 



= -2 N m In (iY m ) - N m - In (N m \) (8) 

m 

where N m and N m are the observed and model differ- 
ential number counts, respectively. Finally, we define 
AC 2 = C 2 — min(C 2 ), w hich will follow the standard x 2 - 
distribution ({Cash 1979) with one degree of freedom (in 
this case the scale height). In Figure |6l we show the ob- 
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TABLE 3 
Ultracool Dwarf Candidates 
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RAT 
r\+V 1 


JJCC 1 
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l,-<600 — J098J 
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/h m s\ 

l ) 
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(mag) 


(mag) 


(mag) 




1 


r\c O Q 0"7 Q 1 

Ub 23 27.31 


f A Q "1 oo n 

— b4 ol zz.U 
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U.o4 ± U.Uo 
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oq a 1 _L n no 
26. 41 ± U.Uz 
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1 O n "1 1 n Ar 

13 Ul 13. U5 


1 oo oo oo o 

+0U 0U 09.0 




1 A C 1 O O /I 

1.45 ± U.z4 
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U.43 ± U.U4 


OO A O 1 O OO 
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MLT 


7 


lb 32 17.1b 


+37 33 32.7 


o oc i n o o 
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O OO 1 O O A 

0.89 ± U.04 


O A O 1 O OO 

U.43 ± U.03 


OO C7 1 O OO 

23.57 ± U.02 


MLT 


8 


09 25 37.82 


+40 01 03.9 


2.77 ± 0.90 




0.86 ± 0.07 


0.19 ± 0.05 


23.7b ± 0.03 


MLT 


9 


08 46 16.63 


+76 53 12.5 




0.72 + 0.34 


1.34 + 0.10 


0.b4±0.04 


24.21 ±0.03 


L 


10 


16 32 21.30 


+37 32 52.1 


2.53 ± 0.77 




0.93 ± 0.08 


0.39 ±0.04 


24.21 ±0.03 


MLT 


11 


15 24 08.81 


+09 55 06.2 




1.84 + 0.70 


0.87 + 0.08 


0.19 ±0.07 


24.b5±0.05 


MLT 


12 


16 31 32.80 


+37 35 53.6 


>3.12 




0.8b ±0.13 


0.19 ±0.10 


24.72 ± 0.06 


MLT 


13 


04 39 21.53 


-53 16 52.0 




0.65 + 0.30 


l.lb ±0.11 


0.49 ± 0.04 


24.73 ± 0.03 


L 


14 


06 23 39.94 


-64 30 58.3 






0.85 ±0.29 


0.35 ±0.19 


24.86 ±0.12 


MLT 


15 


09 14 22.12 


+28 21 34.6 


>3.01 




0.7b ±0.10 


0.18 ±0.07 


24.93 ±0.05 


MLT 


16 


04 39 26.10 


-53 lb 01.1 




>2.92 


0.93 ±0.11 


-0.28 ±0.07 


25.08 ±0.03 


T 


17 


08 19 19.00 


+49 11 02.3 


>2.bl 




0.81 ±0.21 


0.25 ±0.18 


25.28 ±0.11 


MLT 



t Coordinates refer to the J2000 equinox. 

t Spectral types based on equations Jl}, 10, and l3t . 



er ly interpret the star counts. Whil e the dust maps 
of Schlcgel, Finkbc iner. fc Davisl (|1998l) provide the best 
estimate of the total line-of-sight Galactic extinction 
for extragalactic objects, they cannot be directly ap- 
plied to our objects which reside within the Galaxy. 
Instead of parameterizing A(l, b,R), we perform the 
above minimization for both A(£, b, R ) = mag and the 
iSchleeel. Finkbeiner, fc Davisl (|1998t ) value as given m 
Table [TJ As this approach brackets the two extinction 
extremes (the minimum and maximum extinction mod- 
els, respectively), we expect it to indicate the degree to 
which insufficient knowledge of A(£, b, R) is adversely af- 
fecting our results. Therefore we take the average and 
difference between the two extremal dust hypotheses as 
the expected value of the scale height, and its systematic 
uncertainty due to the extinction model, respectively. 
We approximate this as a symmetric uncertainty. How- 
ever we do not expect these objects to be uniformly dis- 
tributed between the two extinction limits. We find that 
the range of scale heights for the minimum and maximum 
extinction models is ±5 pc. 

5.2. Ultracool Dwarf Luminosity Function 

The measured ultracool dwarf luminosity function has 
potentially sizeable uncertainties ~ 30%), largely 

due to the Poisson counting uncertainty of these rare ob- 
jects. To estimate the uncertainty introduced by the in- 
accurate knowledge of the luminosity function, we draw a 
normal random number for each absolute magnitude bin 
with mean and variance from the p ublished results (e.g. 
iCruz et all 120071 iRevle et all 12010) . We recompute the 
vertical scale height for 1000 realizations, and take the 
dispersion of optimal scale heights to be the systematic 
uncertainty. We estimate this uncertainty for our fields 
to be ±30 pc. 

5.3. Malmquist Bias 

For any flux-limited survey, intrinsically bright objects 
are preferentially detected, which biases the mean abso- 
lute magnitude (M) as a function of apparent magnitude 
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Fig. 6. — Ultracool dwarf number counts. The histogram rep- 
resents the number counts from 231.90 arcmin 2 from 39 parallel 
fields and the ERS data in the GOODS-S field. The solid line and 
grey region shows best-fitting model counts for z sc ; = 290 pc and 
the corresponding la uncertainty, respectively. In the upper left, 
we show the AC \z sc {) used to derive these scale height values. 
The peak in the model number counts at J125 ~ 24 mag is not due 
to generic object incompleteness at the faint-end, but rather from 
having integrated completely through the disk. 

served (solid histogram) and model counts (dashed line) 
for the optimal model of z sc i — 290 with the total la 
uncertainty (shaded region — we discuss our systematic 
uncertainties in §[5]), as well as the AC 2 (z sc i) curve for 
A j = mag in the inset (discussed in detail in § 15. ip . 
We compute the random uncertainty on the scale height 
where AC 2 (z sc i) = l to be ±25 pc. 

5. SYSTEMATIC BIASES AND UNCERTAINTIES 

Here we discuss sources of potential systematic uncer- 
tainties that could affect our estimate of the vertical scale 
height. 

5.1. Interstellar Extinction 

As mentioned above, it is necessary to incorporate 
the field- and distance-dependent extinction to prop- 
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(|Malm quist 1936). Under basic assumptions, the correc- 
tion to the absolute magnitudes is given by 

AM = M-M , (9) 



where Mq is the intrinsic absolute magnitude, a is 
the dispersion on the brightnesses from the width of 
the main sequence an d photometric uncertainties (e.g. 
iBochanski et al.H2010T >. and log e = 0.4343. Since our ob- 
served counts (N m ) are very discontinuous due to small 
number statistics, we opt to impose Malmquist bias on 
the model counts (N m ). We adopt a = 0.2 mag which is 
a somewhat more conservative estimate than the typical 
photometric uncertainty of our faintest sources (see Ta- 
ble , which gives shifts on the absolute magnitude of 
— 0.1 < AM < 0.2 magQ. Since this shift is considerably 
smaller than the width of our apparent magnitude bins, 
the bias on the vertical scale height is negligible. 

5.4. Equal-Mass Binaries 

We expect a fraction of our ultracool dwarfs will be 
in binary systems, which if unaccounted for, will tend 
to increase the vertica l scale height measurements (e.g. 
IBochanski etaLl 12010). To assess the properties of po- 
tential binary systems in our sample, we construct a 
grid of simulated images with two point sources placed 
at a range of separations (1 < s < 5 pix and As = 
0.25 pix), total magnitudes (20 < J125 < 25 mag and 
A J125 = 0.25 mag), and Gaussian noise field with mean 
and variance consistent with the parallel fields. For each 
brightness and separation, we generate 1000 realizations 
and catalog the images as described in § 12.31 We find 
that for separations of s < 3 pix SExtractor does not 
detect two distinct point sources, but does recover the 
total flux to ~ 2%. Furthermore, the combined source 
only fails to pass our axis ratio criter i on of (b/a) < 0.7 
for Jus > 23 mag. IBurgasser et al.l (<2007D find that 
most very low-mass stars have physical separations of 
A < 20 AU, which implies that the unresolved bina- 
ries that may escape our cataloging are at a distance 
of 50 < R < 80 pc. Such systems will have an absolute 
magnitude of Mj > 18 AB mag, which corresponds to a 
spectral type t hat is far too cool to have been included 
in our sample ()Hawlev et al.l l2002f) . Therefore we did 
not systematically reject any marginally-resolved bina- 
ries based on our axis ratio criterion. 

Unresolved equal-mass binaries will be 2.51og(2) mag 
brighter than a single star of the same spectral type 
and distance, which will skew the observed counts to 
bright er values and give th e impression of a thinner 
disk (Bochanski et al. 2010). To estimate the magni- 
tude of this effect, we randomly select a fraction of our 
objects (denoted as /&) to be equal-mass binaries. We 
dim these objects by 2.51og(2) mag, duplicate their en- 
tries in the number counts if the dimmed brightness is 
J125 <25.5 mag, and recompute the vertical scale height 
according to § [4] We repeat this procedure 1000 times 
for fb =0.1, 0.2, 0.3, and 0.4 and find that the biases on 
the scale heights are (z^ ary — z sc i) = 5 ± 4 pc, 6 ± 10 pc, 

12 Since our model counts peak at J125 ~ 24 mag, the absolute 
magnitude shifts are not always positive. 



14 ± 12 pc, and 17 ± 12 pc, respectively. The uncertain- 
ties in these biases reflect the distribution of estimated 
scale heights. 

6. DISCUSSION 

With the deep (J 12S < 26 mag) HST/WFC3 parallel 
and pointed fields, we can identify an LO-dwarf out to 
R ~ 3.5 kpc and a TO-dwarf to ~ 700 pc. Since these fields 
are at high Galactic latitudes, they provide constraints 
on the vertical distribution of these intrinsically very 
faint objects. From our compilation of 231.90 arcmin 2 
of HST imaging, we have identified 17 ultracool dwarf 
candidates, whose number counts are consistent with an 
exponential vertical distribution with a scale height of 
z sc i = 290 ± 25 (random) ± 31 (systematic) pc. Our es- 
timate im proves upon p r eviou s results by combining the 
depths of IPirzkal et al.l (120051 ) with the multiple sight- 
lines and area of iRvan et al.l ([2005] ). Additionally, our 
sample likely contains fewer M-dwarf s, owing to the deep 
infrared imaging. For example, the IRvan et al.l (2005) 
work identified dwarfs later than ~M6 from HST/ ACS 
parallels with a single color of (i 1 — z') > 1.3 mag (see 
IBochanski et al.l 120071 for representative SDSS colors of 
M-dwarfs). Since the early M-dwarfs are of consider- 
ably higher luminosity and far more common tha n the 
L-dwarfs, we speculate that the IRvan et a l. ( 2005) sam- 
ple contains many M-dwarfs. Using our derived Galac- 
tic structure model, we estimate that our sample con- 
tains 6.0 ± 2.2, 7.4 ± 2.2, and 2.3 ± 0.3 M8-M9, L, 
and T-dwarfs, re spectively. Had we adopted the weaker 
color criteria of IRvan et ail (|2005t ). we would expect 
27 ± 10 M6-M9 dwarfs. Given these likely differences 
in sample properties, it is not surprising to find possi- 
ble differences in the vertical scale height measurements. 
Finally, we note that the model counts peak around 
J125 — 24 mag, much brighter than the 50% completeness 
limit of J125 — 26 mag (see Table[T]). This peak is not due 
to generic survey incompleteness at the faint-end, but 
rather due to the number density declining faster than 
the volume surveyed. At present, our observed number 
counts do not show or strongly demand such a peak, and 
more dwarfs at J125 ^24 mag are needed to identify this 
critical peak. 

We have used the most recent esti mates of the ul- 
tracool dwarf lum inosity function (e.g. iCruz "etani2"00"7t 
IRevle et al.l |2010T ). which are derived primarily from 
nearby samples (R < 100 pc). Since objects below the 
hydrogen-burning limit are passively-cooling, their bolo- 
metric luminosity strongly depends on their age and ini- 
tial temperature. Therefore, the cooling will introduce a 
non-trivial time dependence on th e luminosity funct ion of 
a population of ultracool dwarfs (Burgasscr 2004c). For 
example, if the majority of these dwarfs are formed at the 
disk midplane and are scattered to these high Galactic 
latitudes by interactions with massive objects in the disk, 
then the luminosity function of these dwarfs is likely dif- 
ferent than the local estimates, particularly if the cooling 
times are shorter than the scattering times. Specifically, 
many of the earliest L-dwarfs will have cooled to become 
later types, resulting in a more "bottom-heavy" luminos- 
ity function with respect to local estimates. 

Additionally, this cooling should tend to make 
early- to mid-L dwarfs a k inematically y ounge r pop - 
ulation than the M-dwarfs iSeifahrt et al.l (e.g. I2010f> . 
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However the high velocity dispersions reported by 
many kinematic studies sugges t ages of 1-6 Gyr 
(e.g. iZapatero Osorio et all l2007t iFahertv et all 120091: 
Seif ahrt et alJ 1201(1) . with a well-established age of ~ 
3 Gyr for the M-dwarfs (e.g. lReid et"al] l2002) . If dwarfs 
immediately below the hydrogen-burning limit are in- 
deed kinematically younger (and have a lower velocity 
dispersion) than the lowest mass main sequence dwarfs, 
then we expect they will be distributed in a thinner 
disk. Yet cooler spectral types will contain a mix- 
ture of old (high mass) dwarfs that have cooled and 
young (low mass) objects. This population will then be 
kinematically older, have a higher velocity dispersion, 
and reside in a thicker disk than the warmest brown 
dwarfs. Therefore we expect to see a gradual deviation 
in the trend of scale height with spectral type, since the 
hydrogen-burning limit does not occur for a fixed spec- 
tral type. With the present data, we find a scale height 
of z sc i = 290 ± 40 pc for a sample of M8-T d warfs, which 
is comparable to estimates fo r mid-M dwarfs (|Juric et al.l 
120081: iBochanski et al.ll2010l ) and is qualitatively consis - 
tent with the kinematic results (e.g. lFahertv et al~ll2009l ). 
However, our estimate is somewhat lower than the ex- 
trapolation of the tr end of scale heig ht with spectral 
type (see Figure 10 of Uuric et al.ll2008r l. More data with 
greatly improved spectral typing is needed to fully con- 
strain these effects. 

Our sample may contain dwarfs as early as ~M8, which 
is a direct consequence of the filter set. These infrared 
colors are determined mostly by a series of molecular 
absorption bands, notably H2O and CH4, which are in 
turn used to define the spectral types (e.g. iKirkpatrickl 



120051 and references therein). Therefore, a cleaner se- 
lection and spectral typing can be achieved by using 
medium or narrow bands whi ch isolate these spectral 
features (e.g. Uones et al.lll994f ). For example, the H2O 
feature at A = 1.34 11m. direc tly maps onto effective tem- 
perature (|Jones et al.lll995l ) and is relatively insensitive 
to surface gravity a nd metallicity (jGorlova et al.l [20031 : 
iWilking et al.ll2004D . Unfortunately the J125- and H 160 - 
bands equally split the H2O feature, diminishing the type 
discrimination of the ( J125 — Hi&o) color. Future surveys 
dedicated to finding ultracool dwarfs could take advan- 
tage of these molecular features for robust identification 
and classification. Fortunately, WFC3 contains a host of 
filters designed to sampl e this H2O absorption feature 
(|Lupie fe Boucarut|[2000T) . specifically F127M, F139M, 
and F153M. Furthermore, the James Webb Space Tele- 
scope (JWST) and its Near-Infrared Camera (NIRCam) 
will be equipped with similar bandpasses, but with a 
significantly larger collecting area facilitating surveys at 
still larger heliocentric distances and search for ultracool 
dwarfs associated with other Galactic components (such 
as thick disk, halo, or bulge). 
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